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The  research  project  "Non  parametric  and  Sequential  Analysis 
of  Life  Testing  and  Reliability  Problems"  started  in  April  1, 

1975  and  was  completed  in  May  31,  1979.  During  this  period  the 
principal  investigator,  along  with  his  collaborators  (who  were 
not  supported  by  AFOSR) ,  considered  a  number  of  problems.  We 
describe  them  below  briefly.  Lists  of  publications  and  techni¬ 
cal  reports  prepared  under  the  project  are  given  in  the  appendix. 


1)  Classification  ^nd  selection  rules  for  the  exponential  population^ 
Let  Hq,  ,  ,Trj^  be  k+1  different  exponential  populations,  j 

The  problem  of  classification  is  to  classify  as  ir^  for  exactly 
oii&  X  (i— 1 , 2  / . . . , k) .  Previously  we  have  developed  classification  /,  ■ - 
procedures  for  the  one  parameter  exponential  distributions.  C’ J 

(Classification  rules  for  the  exponential  populations,  by  A,  P, 

Basu  and  A.  K.  Gupta  (1973)  published  in  Reliability  and  Biometry; 
Statis tical  Analysis  of  life  length ,  ed.  Proschan  and  Serfling,  j 

SIAM  p.  637—650)  •  We  have  extended  many  of  the  results  for  one— 

‘  I 

parameter  exponential  distributions  to  the  case  of  the  two— para¬ 
meter  exponential  distributions  [4j,  Classification  rules  based  i 

the  likelihood  ratio  criterion  for  the  two— parameter  exponential 
populations  have  been  studied.  These  classification  rules  have 
been  extended  to  many  populations  and  to  the  situation  where 
observations  are  censored.  In  many  cases  where  the  likelihood- 
ratio  criterion  does  not  lead  to  appealing  classification  rules, 
ad  hoc  rules  are  proposed.  An  additional  rule  is  considered  from 
testing"  point  of  view,  in  each  case,  the  probability  of 
misclassification  is  derived  exactly  where  the  parameters  are  ‘ 
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known.  The  rules  considered  have  been  shown  to  possess  a  con¬ 
sistency  property.  An  alternative  approach  based  on  Bayesian 
considerations  is  also  explored. 

A  nonparametric  classification  problem  as  described  below, 
has  also  been  considered. 


2)  Nonparametric  analysis  of  some  reliability  problems. 

Suppose  a  random  sample  units  with  failure 

time  distribution  F(x)  is  put  on  test.  Similarly,  let  a  second 
sample  **  units  from  G(y)  be  tested.  To  test  the 

hypothesis 


Hq:  P  =  G 

against  the  alternative  that  they  are  different,  based  on  the  first 
r  ordered  observations  (failure  times)  out  of  a  combined  sample 
size  m+n,  the  following  statistic  is  considered. 


^  1 

£  (nm.-mn. )  +  irCN-r-l)  (nm  -mn  )  (r^N) 
1  1  2  r  r 


where  m^  and  n^  are  the  nxsnber  of  x  and  y  failures,  respectively, 
among  the  first  i  ordered  observations  of  the  combined  (x,y)  sample. 

Properties  of  and  its  various  k-samp..^  extensions  are  re¬ 
viewed.  Suitedale  modification  of  is  used  to  classify  an  unknown 
population  into  one  of  two  given  populations.  The  classification 
rule  is  shown  to  be  consistent  and  its  asymptotic  relative  effi¬ 
ciency  with  respect  to  a  classification  rule  for  the  normal  dis¬ 
tribution  is  computed.  The  problem  of  estimating  the  reliability 
fionction  P(Y<X)  is  also  discussed. 
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3i ^ paves ian  test  for  increasing  failure  rate^ 

In  situations  where  a  statistical  analysis  is  carried  out  on 
life-testing  data,  it  is  often  useful  to  study  the  failure  rate 
of  the  population.  The  failure  rate,  q(t,x),  of  a  population 
over  the  time  interval  (t,t+x)  is  defined  as  the  ratio  of  the 
probability  that  an  item  selected  at  random  from  the  population 
will  fail  in  the  time  interval  (t,t+x) ,  given  that  it  stirvived 
up  to  time  t,  divided  by  x.  It  has  been  found,  both  by  theoretical 
and  practical  investigators,  that  many  of  the  probability  distri¬ 
butions  useful  in  life-testing  have  either  a  monotone  increasing 
or  monotone  decreasing  (or  constant)  failure  rate.  That  is, 
q(t,x)  is  a  monotone  function  of  t  for  fixed  x.  Several  tests 
^have  been  developed  by  others  to  test  if  a  population  has  mono¬ 
tone  failure  rate.  In  this  article  inference  for  q(t,x)  is  de¬ 
veloped  from  a  Bayesian  viewpoint.  The  failure  rate  for  a  future 
observation  is  obtained  based  on  an  existing  data  set.  Bounds 
on  the  posterior  probabilities  of  a  distribution  with  increasing 
or  decreasing  failure  rate  are  also  obtained.  The  methods  can 
be  used  with  random  sauries  or  multiply  censored  samples. 

4)  Estimation  of  the  restricted  scale  parameter  of  the  exponential 
distribution,  ^ 

^l'^'****^n  of  n  items  put  on  test  and 

the  density  function  of  the  underlying  life  distribution  be 

I 

given  by  “  ’  '  _ _ ’ 

f(x,e)  *  I  e’""/®  ,  x>0  . 

Suppose,  based  on  past  data,  it  is  known  that  8>0  ,  where  8  is 

o  o 

known.  Such  a  situation  arises  naturally  in  mauiy  physical  problems, 
where  it  may  not  be  possible  to  make  specific  statement  about  in- 
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dividual  observations,  but  usually  it  is  possible  to  guarantee  a 
mininuim  value  for  the  mean  life.  Thus  in  a  life  testing  context 
there  might  be  an  early  failure  due  to  production  error,  or  due  to 
the  presence  of  some  outliers,  but  overall  mean  life  may  be  known 
to  be  greater  than  some  minimum  value. 

Several  estimators  of  0  are  being  considered  and  their  prop¬ 
erties  are  being  studied. 

5)  Nonparametric  tests  for  independence . - ' 

The  object  of  this  study  is  to  consider  rank  tests  for 
independence  with  possible  applications  in  life  testing’ and  bio¬ 
medical  problems.  The  work  has  been  carried  out  in  two  different 
directions. 

The  first  direction  is  to  investigate  properties  of  rank 
tests  under  appropriate  models.  Under  alternatives  of  bivariate 
normality,  performance  of  snch  tests,  specifically  of  Kendall’s 
T,  Spearman's  p,  and  the  normal  scores  statistic,  is  well  kno^m. 

In  the  area  of  life  testing  more  interest  is  given  to  distri¬ 
butions  of  the  bivariate  exponential  type.  Of  these,  those 
possessing  the  lack  of  memory  property  are  of  special  interest, 
for  example,  the  models  of  Marshall  and  Olkin  (BVE)  and  Block 
and  Basu  (ACBVE). 

Testing  for  independence  under  the  BVE  distribution  is 
shown  to  reduce  to  that  of  the  ACBVE.  For  the  ACBVE  model. 

Pitman  asymptotic  relative  efficiencies  are  calculated  for 
Kendall's  r,  Spearman's  p,  the  maximum  likelihood  statistic, 
and  the  statistic  providing  the  uniformly  most  powerful  test. 

An  appropriate  exponential  scores  statistic  and  the  locally  most 
powerful  rank  statistic  are  derived.  Approximate  large  sample 
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power  is  calculated,  and  small  sample  power  is  investigated 
using  Monte  Carlo  methods. 

Similar  but  less  extensive  work  is  done  for  a  special  case 
of  the  tri.variate  extensions  of  the  BVE  and  ACBVE  distributions. 
Two 'parametric  tests,  one  the  likelihood  ratio,  are  derived  as 
are  an  exponential  scores  statistic  and  the  locally  most  power¬ 
ful  rank  statistic.  Power  studies  are  carried  out  for  these 
statistics  and  for  a  tri variate  generalization  of  Spearman's  p. 

The  general  trivariate  extensions  with  the  restriction  of 
identical  marginals  are  discussed  along  with  the  resulting 
likelihood  ratio  test.  The  results  are  easily  extended  to  the 
k-variate  situation. 

A  second  direction  has  been  to  use  suitably  modified  rank 
tests  to  analyze  censored  data  arising  from  life  testing  and 
other  biomedical  situations.  An  extensive  computer  study  is 
done  to  investigate  the  behavior  of  several  statistics  which 
are  modifications  of  Kendall's  t  adjusted  for  censoring.  Power 
simulations  are  done  for  both  the  ACBVE  and  bivariate  normal 
distributions  subjected  to  random  censoring  in  both  variables. 
These  statistics  are  then  applied  to  test  for  independence  in 
heart  pacemaker  patient  survival  data. 

-^6)  Estimate  of  reliability  in  the  stress-strength  model^  ' 

Sixppose  Y  is  the  strength  of  a'  component  which  is  subject 
to  a  stress  S.  Then  the  component  fails  whenever  X  i  Y,  and 
there  is  no  failure  when  X  <  Y.  The  problem  of  estimating  the 
reliability  function 
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R  =  P(X  <  Y) 

is  considered.  A  survey  of  available  results  is  presented  and 
some  new  results  are  considered  when  X  and  Y  follow  gamma, 
exponential  or  Weibull  distribution.  The  case  when  X  and  Y 
follow  a  bivariate  exponential  distribution  is  also  considered. 

7)  Estimates  of  reliability  for  k-out-of-m  systems. 

■  Let  S  be  a  system  of  m  independent  components  s^,S2» 
...,Sjj^  such  that  the  system  operates  if  and  only  if  at  least 
k  of  these  m  components  successfully  operate.  Component 
lives  are  assumed  to  follow  independent  exponential  distributions 
with  tinknown  scale  parameters.  The  minimum  variance  unbiased 
estimator  of  the  reliability  of  the  above  system  at  a  given 
mission  time  t  is  obtained.  The  performance  of  the  estimator, 
when  the  component  lives  ares  identically  distributed,  is  compared 
with  the  corresponding  maximum  likelihood  estimator  for  both 
large  and  small  samples.  Finally,  a  confidence  interval  for 
system  reliability  is  obtained. 

8)  Identif iability  of  the  multinormal  and  other  distributions 

under  competing  risks  model"!  ^ 

Let  Xj^ ,X2f . . .  »Xp  be  p  random  variables  with  joint  distribu¬ 
tion  function  FCXj^, . . .  ,Xp)  .  Let  Z  -  min(Xj^,X2» . . .  ,Xp)  and 
^  i  if  Z  “  X^-.  In  this  paper  the  problem  of  identifying  the 
distribution  function  F(Xj^, . . . ,Xp) ,  given  the  distribution  Z 
or  that  of  the  identified  minimum  (Z,I),  has  been  considered 
when  F  is  a  multivaraiate  normal  distribution.  For  the  case 
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p  ■  2,  the  problem  is  completely  solved.  If  p  *  3  and  the  dis¬ 
tribution  of  (Z,I)  is  given,  we  get  a  partial  solution  allowing 
us  to  identify  the  independent  case.  These  results  seem  to  be 
highly  nontrival  and  depend  upon  Liouville's  result  that  the 
(Univariate)  normal  distribution  function  is  a  nonelementary 
function.  Some  other  examples  are  given  including  the  bivariate 
exponential  distribution  of  Marshall  and  Olkin,  Gumbel,  and  the 
absolutely  continuous  bivariate  exponential  extension  of  Block 
and  Basu. 

Concluding  Remark 

The  present  project  solved  a  number  of  problems.  It  also 
raised  new  ones.  Some  of  these  will  be  studied  further  by  the 
principal  investigator  in  the  future. 
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